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COMMON ENVELOPE: ENTHALPY CONSIDERATION 

N. IVANOVA 1 , S. CHAICHENETS 2 
Draft version January 14, 2013 

ABSTRACT 

In this Letter we discuss a modification to the criterion for the common envelope (CE) event to result in envelope dispersion. 
We emphasize that the current energy criterion for the CE phase is not sufficient for an instability of the CE, nor for an ejection. 
However, in some cases, stellar envelopes undergo stationary mass outflows, which are likely to occur during the slow spiral-in 
stage of the CE event. We propose the condition for such outflows, in a manner similar to the currently standard acEA-prescription 
but with an addition of P/p term in the energy balance equation, accounting therefore for the enthalpy of the envelope rather than 
merely the gas internal energy. This produces a significant correction, which might help to dispense with an unphysically high 
value of energy efficiency parameter during CE phase, currently required in the binary population synthesis studies to make the 
production of low-mass X-ray binaries (LMXBs) with a black hole companion to match the observations. 
Subject headings: binaries: close — stars: evolution — X-rays: binaries 



1. INTRODUCTION 

The most important event in the formation and evolu- 
tion of most close interacting binaries is the so-called com- 
mon envelope (CE) event, during which the components of 
a binary system are engulfed by a common gaseous enve- 
lope, and the resulting interaction dramatica l ly shrinks their 
orbit (iBisnovatyi-Kogan & Syunyaevl 11971b lOstriker! 119761: 
lPaczynskilll976h . Depending on the envelope structure and 
companion masses, the envelope is ejected to leave behind a 
close binary, or the two stars merge. Dividing the parame- 
ter space into binaries that survive CE and those that do not, 
and determining the final separation for the former, is criti- 
cal for calculating formation rates of low mass X-ray binaries 
(LMXBs), 7-ray bursts, as well as for LISA and LIGO events 
(Belczynski et alJl2008h . 

In the standard treatment of CE outcomes, the final sep- 
aration of the b i nary is determined v ia "energy formalism" 
dWebbinklll984t iLivio & Sokerj|1988l) . in which the binding 
energy of the (shunned) envelope is equated to the decrease in 
the orbital energy E or \,: 
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Here a\ and Of are the initial and final binary separations, m\ 
and mi are the initial star masses and wi iC is the final mass of 
the star that lost its envelope. 

£bind is assumed to be the energy expense needed to remove 
the envelope to infinity and is commonly adopted to be the 
sum of the potential energy of the envelope and its internal 
energy. To characterise the donor envelope central concen- 
tration and simplify calculations, specifically for population 
synthesis, a parameter A was introduced: 
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Here mi |e is the mass of the removed giant envelope and R[ is 
the radius of the giant star at the onset of CE, e is the specific 
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internal energy and $>(m) = -Gm/r is the (gravitational) force 
potential (or specific potential energy). isA.bind can be found 
directly from stellar structure for any accepted core mass. 

Another parameter, acE, is introduced as a measure of the 
energy transfer efficiency from the orbital energy into enve- 
lope expansion, and the balance of energy is written as 
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Many authors choose to accept ac E A ~ 1> since for many (at 
least low-mass) stars A = 1, as can be found from detailed 
stellar structure, and chce is bounded above by 1 . 

The simplicity of the standard prescription resulted in its 
popular use in the binary population synthesis calculations. 
However, once accurate values of A from stellar structure cal- 
culations were determined, this approach has shown incon- 
sistencies with the observations, especially large for the for- 
mation of black hole LM XBs: in massive giants A <C 0.1 
dPodsiadlowski et aT]|2003l) . and it has been shown that with 
cuce < 1 only an intermediate-mass companion could avoid 
a merger; this challenges the formation of a lo w-mass X-ray 
binary with a low-mass companion in general (Just ham et al.l 
2006). 

Here, we revise the energy requirements necessary to force 
a common envelope to disperse and discuss its possible appli- 
cation for LMXBs formation. 

2. ENTHALPY CONSIDERATIONS 
2.1. Total energy and instability 

While the defini tion of the gravitation al binding energy is 
unequivocal, (e.g., Chandrasekhar 1939), there exists no au- 
thoritative source defining what expression should be used as 
the binding energy of a gaseous sphere in the sense of § 1. 
To clarify, in treating the common envelope problem we are 
interested in the additional energy to be deposited in the en- 
velope in order to disperse material to infinity. The 'industry 
standard' at the moment is to use the total of the gravitational 
U and the gas thermal E energies, which, as we shall argue in 
this section, is not correct. 

As an illustrative limiting case, consider a star with the 
positive total energjQ Wo|f=o > to begin with, which 

3 with zero of energy defined for all material evacuated to infinity, and 
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is "kinetically stable"; i.e. an energy barrier has to 
be overcome between the bound a nd unbound states 
dBisnovatvi-Kogan & ZerDovicb] [i967). It is clear that in 
such cases, the additional energy A required to unbind the 
star is the magnitude of the energy barrier, not the (unphysi- 
cal)-W <0. 

The secular stability of a star against small adiabatic pertur- 
bations is not defined by the sign of the total energy W, but 
rather by the variational conditions: an equilibrium state is an 
extremum of the total energy (<5W) ac j = 0; and secularly stable 
configuratio ns are then found at local minima: (<5 2 W) a d > 0. 
dChiuli 19681) . 

If the first adiabatic exponent T[ is approximately a con- 
stant throughout the star, the secular stability criterion reduces 
to the condition dChiull 19681) 



(4) 



Here P is the pressure and p is the density. 

For a one-zone model of a stellar envelope, a li nearized ver- 
sion of the above condition, the Baker's model ( Baker! [T9661) 
gives a similar criterion for the volume-averaged adiabatic ex- 
ponent Tj > 4/3 in the envelope 

On the other hand, the virial theorem, applied to the en- 
tire star, gives us the condition for the total energy W < 
in the case of cons tant third adiabatic exponent T^, as (e.g., 
lHansen et al.ll2004T) : 
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Here T is is the temperature. We note that this simple form 
depends crucially on Newton's third law applied to every pair 
of particles inside the star, and thus needs to be substantially 
modified if we only consider the envelope. 

We reiterate the importance of distinguishing these condi- 
tions: instability does not have to occur in a state with positive 
total energy, and Ti does not have to coincide with 1^; they 
are related as 

In particular, the ionization zones, where V a d = (§|^) ad < 
0.4 and can become as low as 0.1, have local value of T$ < 
4/3 while Ti > 4/3. 

It is unfortunately non-trivial to find the volume-averaged 
T[ in the giant envelope after depositing some heat in it. How- 
ever, it is clear that the criterion of the total energy in the en- 
velope to be W = is not applicable for the envelope stability, 
nor is it a sufficient condition for the envelope to be dispersed 
to infinity. We will consider instead the stability of a stellar 
envelope towards creating outflows. 

2.2. The condition for the envelope to outflow 

Since we are interested in the energy requirements, it is nat- 
ural to consider a lower bound on A: all orbital energy is con- 
verted into heat, and the velocity of ejecta at infinity is zero. 
How precise these bounds are is not of great concern here and 
is a separate problem on its own, although we expect them 
to be reasonably close: the time scale for viscous friction be- 
tween differentially rotating regions can exceed the spiral-in 

pressure set to zero 



time, but most angular momentum is lost in the outer enve- 
lope while most of the energy is released at closer orbits, so 
the overall effect of angular momentum conservation can at 
least be constructed to be small. 

As any other thermodynamical system undergoing a steady 
process, the material in a star obeys the first law of thermo- 
dynamics, which states that the change of the internal energy 
comes not only from the heat transferred into system, but also 
from the work done by the system: 



SE = 5Q-P8 



(7) 



Applying t his equation for each m a ss shell in the envelope , 
(see, e.g., IColgate & White! 1 19661: iKutter & Sparks] 119721) . 
gives the energy conservation equation in Lagrangian coor- 
dinates as 



d (u 1 T \ 8(1/ p) n 
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Here u is the velocity. 

The initial condition is that at the start the mass shells in the 
envelope are not moving, but every shell in the envelope re- 
ceived some heat, in this case deposited from the companion's 
orbital contraction. The heat has some arbitrary distribution 
Sq(m) > in the envelope (Sq(m) is per mass unit), such that 
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{m)dm = Q. 
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A given lagrangian shell, once it has been heated 
and started expansion, will have reached the point 
of no return in its expansion when it s E > (e.g., 
i Bisnovatvi-Kogan & Zel'DovicrJ 119671 : ISparks & Kuttel 
11971 

PI p 

((fy(m)+#+e+-) s tart = (x" 2 + *+e+-)ex P = S = const (10) 
p 2 p 

In this form, the equation can be also recognized as version of 
the Bernoulli equation. 

The quantity E, for stability analysis purposes, is more 
important than the total energy: bipolitropic stellar model 
that has E < in its envelope but positive total energy is 
metastable as a whole, even though the run-away to in- 
finity is energetically allowable; but a star with E > 
in its envelope will be always quasi- s teadily outflowing 
(IBisnovatvi -Kogan & Zel'Dovichl 119671 : iBisnovatyi-Kogarl 
2002). It also has been noted that it is a general feature of low- 
mass giants during double shell burning to establish E > in 
larger and large parts of their envelopes as they evolve (with 
He shell providing large energy inflow to the envelope with 
each He-shell flash), and could possibly being responsible for 
envelop e outflows and/or ejecti ons with planetary nebula for- 
mation (ISparks &Kutterlll972l) . 

In line with res ults of past stability anal yses of when a star 
starts an outflow dBisnovatvi-K ogan 2002, and the references 
therein), we suggest that once a part of the common envelope 
has obtained positive E, it will start outflowing, notwithstand- 
ing the sign of the envelope's total energy. A stronger crite- 
rion would be to require the entire envelope to have positive 
E. Under assumption of the minimum energy requirement 
(velocity of gas at infinity is zero), we can write the energy 
conservation for a whole envelope as 
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FIG. 1. — \t and h inside 2 Mq giant. The "boiling pot" zone has excess 
energy. Energy is the specific energy in erg per g. 

/■ surface / P(m)\ 
Q+ *(m) + e(m) + -— - )dm = 0. (11) 

./core V P(' n ) J 

The timescale for start such quasi-steady surface outflows 
is one on which the envelope redistributes the dumped heat, 
i.e. the thermal timescale of the envelope - recall our earlier 
note about most energy being deposited in lower orbits. This 
time is about few hundre d years (e.g, it is about 1000 years 
for a 2QMq, see Fig.l in Ilvanoval l201 1). This makes the ap- 
pearance of mass outflows natural during the self-regulating 
spiral- in stage; such a sta ge could last for up to a thousand 
years (Podsiadlowski 2001); but they likely will not take place 
if the common envelope event occurs on the dynamical time, 
as is the case e.g., in the case of a physical collision, or if the 
swallowed companion is too small to expand the giant enve- 
lope and establish a self -regulated slow spiral-in phase. 

Note, that 

• due to the presence of a non-negative term P/p, this 
condition occurs before the envelope's total energy be- 
come positive; 

• the master equation (0 of the ' A-formalism' is a version 
of eq.dlOl) where all the velocities, as well as the work 
P/p are simply neglected. 

The term P/p is of the order of magnitude of e; for exam- 
ple, for an ideal gas, with no radiation pressure and ionization 
taken into account, Pj p = 2/3 e. The quantity h = e + P/p is 
generally known as enthalpy. 

In line with the classical 'A-formalism' and for ease of com- 
parison, we introduce: 

/•surface /~> 

E hMnd = - mm)+h(m))dm= ' u (12) 

Jcore Ah-Kl 

Following ISparks & Kutterl d 19721) in considering the de- 
tailed stellar models of giants with respect to their E, we 



also note that normal single giants have a "boiling pot" zone 
(BPZ), where E(;«) = ^(m) + h(m) > due to high value of 
P/p (Fig. [T]). Once sufficient amount of the "lid" - the star 
matter above this zone is removed, all material with positive 
E could freely stream away without the need to convert any 
additional mechanical energy. The outer boundary of the BPZ 
almost does not change during the giant stage, and the inner 
boundary does not change once the convective envelope is es- 
tablished (e.g., in case of low-mass red giants, RGs). The 
bottom of the BPZ roughly (but not always exactly) coincides 
with the bottom of the outer convective zone. Illustratively, 
the mass contained in this zone is: ~ 0.5M Q for 2M© RG, up 
to 7.7M Q for a 30M Q giant evolved without wind mass loss, 
and about 3 times less for a 3OM giant evolved with mass 
loss. We anticipate that the presence of the boiling pot zone 
could be a key to understand why RGs are expanding when 
losing their mass. 

3. COMPARISON OF E a ,b IND AND E H ,BIND IN STELLAR MODELS 

Just like A, Ah can be found from detailed stellar calcula- 
tions. For the purpose of this comparison, we adopt t he def - 
inition of the core mass for both cases as in llvanoval (1201 ll) : 
the post-CE core is the point with maximum compression P/p 
inside the Hydrogen burning shell, mi iCp . Such a core will not 
re-expand after the envelope loss, and it was shown that in 
most of the cases it is energetically beneficial to remove the 
envelope down to exactly mi. cp . 

In Table 1 we show the results for giants of differ ent masses . 
We used the set of stellar models described in llvanoval (120 1 lh . 
We find that in most giants the ratios of Ah/A are from ~ 2 
to ~ 5, where the largest ratios, for a giant of the same mass, 
are for an earlier giant. The final separations allowed by the 
enthalpy-consideration are then larger by the same few times. 

The minimal companion mass - in the sense that no pos- 
sible stripped core value exists such that a binary would not 
merge - changes by several times as well (Table 1). Specifi- 
cally, we want to point out the difference between the enthalpy 
consideration and standard prescription for massive giants: 
while acEA-formalism would predict the minimum surviving 
companion mas s of about several M w , in line with the prob- 
lem raised by Justh am et al.1 ([2006), the enthalpy-formalism 
allows for an low-mass companion to survive and form in a 
future a black hole LMXB. 

4. CONCLUSIONS 

In this Letter we considered the termination of a common 
envelope event at the moment it establishes a quasi-stationary 
mass outflow that would reach a point of no return. Such out- 
flow develops during the slow self-regulating spiral-in phase 
of the common envelope event. We showed that if the com- 
mon envelope is escaping the binary in this way, neglecting 
the P/p term in the standard «ceA energy conservation pre- 
scription is too gross. If the enthalpy rather than internal en- 
ergy is calculated in the energy balance, it makes a crucial 
difference in giants. 

While we anticipate that the estimate of a proper energy re- 
quirement is paramount for the formation rates of all kinds 
of post-CE binaries, we especially emphasize the importance 
of this effect in the case of CE with a massive giant. When 
the standard energy formalism predicts that only a several so- 
lar mass donor could survive CE, and so the formation of a 
low-mass X-ray binary with a black hole accretor is forbidden 
unless acE exceeds 1, the enthalpy-consideration naturally al- 
lows for a low-mass companion survival. 
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TABLE 1 

CE OUTCOMES COMPARISON. 
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A post-CE mass is adopted to be the divergence point mi >cp , for a comparison 
is shown the mass of the hydrogen-exhausted core mi x (where X < 10" 10 ). 
A and A/, connect the energies required to eject the envelope with their 
parametrizations in two formulations (eqs. 2 and 12). mi \ and mif, are the 
minimum companions' masses that could survive a CE event, where mi \ is 
using standard prescription and wi2.f1 is with enthalpy-consideration, mi are 
the current donor masses and m zams are donor masses at the zero-age main 
sequence (if different from mi). Si are the current donor radii. All masses 
are in Mq , R[ is in Rq . 
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